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Determining relationships between different types of quantum correlations in open composite quantum sys¬ 
tems is important since it enables the exploitation of a type by knowing the amount of another type. We here 
review, by giving a formal demonstration, a closed formula of the Bell function, witnessing nonlocality, as a 
function of the concurrence, quantifying entanglement, valid for a system of two noninteracting qubits initially 
prepared in extended Wemer-like states undergoing any local pure-dephasing evolution. This formula allows 
for finding nonlocality thresholds for the concurrence depending only on the purity of the initial state. We then 
utilize these thresholds in a paradigmatic system where the two qubits are locally affected by a quantum envi¬ 
ronment with an Ohmic class spectrum. We show that steady entanglement can be achieved and provide the 
lower bound of initial state purity such that this stationary entanglement is above the nonlocality threshold thus 
guaranteeing the maintenance of nonlocal correlations. 


I. INTRODUCTION 

Quantum nonlocality and entanglement identify two kinds 
of quantum correlations that are at the foundations of quantum 
mechanics and play an essential role in quantum information 
theory HM- Although for pure states of a bipartite system 
the two concepts coincide, for mixed states the presence of 
entanglement does not imply in general nonlocality, the lat¬ 
ter being a stronger quantumness property than entanglement 
Eia. In many scenarios a quantum system exists in mixed 
states, particularly during its dynamics due to the interaction 
with the surrounding environment which induces decoherence 
and therefore mixedness mM- Moreover, the environmen¬ 
tal effects are usually detrimental for entanglement and non¬ 
locality in a configuration of separated noninteracting qubits 
embedded in independent local environments 100. Such a 
situation is typical for quantum communication and informa¬ 
tion protocols where the supports of the single-particle wave 
functions are centered around sufficiently far locations or the 
particles are separated by an energy barrier large enough, as 
in experiments with photons in different optical modes or with 
strongly repelling trapped ions II9HI3. It is thus useful to find 
dynamical equations that relate these two kinds of quantum 
correlations, for instance in order to acquire information on 
the presence of nonlocality just by looking at the quantum en¬ 
tanglement of the overall system. 

The possible existence of closed relations between quanti¬ 
fiers of entanglement and nonlocality is a subject of special 
interest in dynamical contexts EHia. Entanglement can be 
quantified for an arbitrary two-qubit state by the concurrence 
C{t) iflTl . Nonlocality is instead unambiguously identified 
if a Bell inequality is violated. The Bell function B, as de¬ 
fined by the Clauser-Horne-Shimony-Holt (CHSH) form of 
the Bell inequality 10, can be then used to seek whether the 
system exhibits nonlocal correlations which occur with cer¬ 
tainty if > 2. An important goal within this context is to 
prove the existence of threshold values of concurrence which 
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ensure nonlocal quantum correlations during a given evolu¬ 
tions whenever C (f) stays above that threshold. In fact, apart 
its basic interest, this behavior would permit to utilize the sys¬ 
tem state for quantum information processes relying on nonlo¬ 
cality, as device-independent and security-proof quantum key 
distribution protocols 0111 . It has been recently reported 
that Bell function and concurrence satisfy a closed dynami¬ 
cal relation in the case of general pure-dephasing evolution 
for each qubit, provided that the qubits are initially prepared 
in the class of extended Werner-like states m. This result 
has been then utilized to verify the preservation of nonlocality 
by dynamical decoupling pulse sequences in the presence of 
low-frequency noise typical of the solid state ifT^ . 

The aim of this work is to review the pure-dephasing rela¬ 
tion between concurrence and Bell function, giving a formal 
demonstration of it. Successively, we apply the result to a 
paradigmatic system made of two independent qubits locally 
interacting with pure-dephasing noise sources having Ohmic 
class spectrum. By exploiting the trapping of the single-qubit 
coherence which occurs in such a system ll 20 l . we investigate 
the possibility to achieve a stationary entanglement which en¬ 
sures the presence of nonlocal correlations. In particular, we 
find the lower bound of the purity of the initial two-qubit state 
such as to guarantee that the stationary concurrence is above 
the nonlocality threshold. 


II. PRELIMINARIES: BELL EUNCTION AND 
CONCURRENCE 

In order to witness nonlocality we use here the CHSH 
form of Bell inequality, which is the most suitable for an 
experimental test of nonlocality in bipartite quantum sys¬ 
tems composed by two-level (spin-like) systems, or qubits 
0 . Let the operator Os = Os{ds,(l>s) be a spin ob¬ 
servable with eigenvalues ±1 associated to the spin parti¬ 
cle S = A, B defined by Os = O 5 • as, where O 5 = 
(sin 05 cos (/)s, sin 05 sincos 0 s) is the unit vector indi¬ 
cating an arbitrary direction in the spin space and as = 
(erf, ay, erf) is the Pauli matrices vector. The CHSH-Bell in¬ 
equality associated to a two-qubit state p is B{p) < 2, where 
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B{p) is the Bell function defined as ll^ l2TI 

B{p) = \{OaOb) + {OaO'b) + {O'aOb) - (OaOb)I, (D 

where (OaOb) = is the correlation function 

of observables Oa, Ob and O'g = Os{d'g, (p’g). If, given the 
state p, it is possible to find a set of angles { 0 ^, 0 ^, 0 ^, 
and such that the CHSH-Bell inequality is 

violated (B{p) > 2), then the correlations are nonlocal or, 
in other words, the system state exhibits nonlocality. A pro¬ 
cedure to obtain the maximum of ,B(p) and the corresponding 
angles for an arbitrary two-spin- 1/2 mixed state is well known 
l22l . Using it, one finds B^naxip) = 2-y/maXj>fc{Mj -f Uk}, 
where j, k = 1, 2, 3 and Uj’s are three quantities depending 
on the state. 

Entanglement for an arbitrary state p of two qubits is quan¬ 
tified by concurrence Cl El 

Cab = C{p) = max{0, a/Ai" - a/A^ - \/^ - (2) 

where Xi (i = 1,..., 4) are the eigenvalues in decreasing or¬ 
der of the matrix p{ay 0 <Jy)p*{ay 0 ay), with ay denoting 
the second Pauli matrix and p* corresponding to the complex 
conjugate of the two-qubit density matrix p in the standard 
computational basis S = {|1) = |11),|2) = |10), |3) = 
|01),|4) = |00)}. 

Let us take now the two-qubit states whose density matrix, 
in the basis B, has a X structure as 

( Pit 0 0 pi4 \ 

0 P 22 P23 0 I 

0 Ph ft, 0 ■ 

p *4 0 0 P44 / 

This class of states is sufficiently general to include Bell states 
(pure two-qubit maximally entangled states) and Werner 
states (mixture of Bell states with white noise) ll23l . A re¬ 
markable aspect of the X states is that, under various kinds 
of dynamics, the initial X structure is maintained during the 
evolution ll23l l24l . Using the above criterion to obtain ,Bniax 
Ea, the three quantities it’s, in terms of the density matrix 
elements, are found to be 

Ui = 4(|pi 4| -f |p23|)^, W2 = (pil + P44 — P22 — P33)^) 

^3 = 4(|pi4| — |P23|)^5 (4) 

as already reported 1251 . Being ui always larger than u^, 
the maximum of Bell function for X states results to be 
Buiax (p) = 2 y^ + maxj= 2 , 3 {'Wj}, that is 

^max = max{ 61 ,^ 2 }, Bi = 2y'ui -f U 2 , 

B 2 = 2\/ui+U3. (5) 

The expression of concurrence for X states is given by ||71 

Cp{t) = 2max{0, |p23(t)| - \/pii{t)p44{t), 

\Pli{t)\ - \/p22{t)p33{t)}- (6) 


III. RELATION BETWEEN BELL FUNCTION AND 
CONCURRENCE UNDER PURE-DEPHASING 

In this Section we review the relation between quantifiers of 
entanglement and nonlocality for two noninteracting qubits, A 
and B, locally subject to any pure-dephasing interaction with 
the environment. We start by defining the general form of 
the Hamiltonian and the relevant initial states of the two-qubit 
system. 

The dynamics of each qubit is governed by the pure- 
dephasing Hamiltonian (H = 1) 

H = —(Ti; + ^az + Hr, (7) 

where U is the Bohr frequency of the qubit and x represents 
an arbitrary environmental operator coupled to the same qubit. 
The free evolution of the environment is included in Hr. The 
overall Hamiltonian of the two-qubit system is thus iTtot = 
Ha + Hr, with each Hr (S = A, B) given by Eq. 0 . 

The two qubits are supposed to be prepared in an EWL state 

HH 

Pi =r |la) (lal + ^-^ 14 , P 2 =r \2a) {2a\ + ^—^t4, 

( 8 ) 

where the pure parts |la) = a| 01 ) -f 6 | 10 ) and \2a) = 
a 100 ) + & 111 ) are, respectively, the one-excitation and two- 
excitation Bell-like states with -f = 1 (a is real) 
and I 4 is the 4x4 identity matrix of the two-qubit Hilbert 

g ice. The EWL states are a subclass of the X states of Eq. 

and when a = b = l/v/2 they reduce to the Werner 
states, a subclass of Bell-diagonal states El. The purity pa¬ 
rameter r G [0, 1] is a measure of the purity of EWL states 
which is given by P = Tr(p^) = (1 -f 3r^)/4. The two 
EWL states of Eq. 0 have the same value of concurrence 
CpA^) = C'p 2 ( 0 ) = 2max{0, (|a&| -f l/4)r — 1/4}. Initial 
states are thus entangled for r > r = (1 -f 4|a6|)“^. 

We can now claim the following theorem. 

Theorem 1. For noninteracting qubits starting from an en¬ 
tangled EWL state (r > f) and locally subject to any pure- 
dephasing noise, the Bell function and the concurrence satisfy 
the dynamical relation sms 

B{t) = 2\/r2 + [C{t) + (1 - r)/2]2. (9) 

Proof Since the two qubits are noninteracting, the evolu¬ 
tion of entanglement and nonlocality can be simply obtained 
from the knowledge of single-qubit dynamics ll24l . Under 
a pure-dephasing evolution, for each qubit the diagonal el¬ 
ements of the density matrix in the eigenstate basis remain 
unchanged. The single-qubit coherences evolve in time as 
qs{t) = Poi(^)//°oi( 0 )’ Ihe explicit time dependence being 
specified by the environmental properties and the interaction 
term. If the system is subject to pure-dephasing only, the X 
form of the density matrix is maintained during the evolu¬ 
tion. In particular, diagonal elements remain constant while 
off-diagonal elements evolve in time. They are related to the 
single qubit coherences by ||2l|24l p 23 {t) = P 23 {^)qA{t)q*B{t) 
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for the initial state pi and pii{t) = pi 4 { 0 )qA{t)qB{t) 
for p 2 - From Eq. •mi- we immediately obtain the concur¬ 
rences at time t for the two initial states given in Eq. 0 as 
Cp,(f) = 2 max{ 0 , \p23{t)\ - Vpii( 0 )p 44 ( 0 )} andCp^{t) = 
2max{0, \pi 4 :(t) \ - a/p22 (0)/333(0)}. Eor the pure-dephasing 
evolution, we then have Cp^ (t) = Cp^ (t) = C (f) with 

C{t) = 2inax{0,ra\/l^^aP\qA{t)qBit)\-{l-r)/4:}. (10) 


We now turn to nonlocality. Eor independent qubits subject to 
local pure-dephasing noise, the two functions Bi, B 2 defined 
in Eq. Q have the same form for the initial EWE states of 
Eq. (|^ and are given by 

Bi{t) = 2\/r2 -f 4r2a2(l - a‘^)\qA{t)qBitW, 

B2(t) = 4:^/2r\a\VT^^\qA{t)qB{t)\. ( 11 ) 


Since a,r,\qA{t)\,\qB{t)\ S [0,1], it is easily seen that 
Bi{t) > for any t: the maximum Bell function 

is thus B{t) = Bi{t). We now observe that nonzero 
entanglement (C(f) > 0 ) is necessary in order that the 
state exhibits nonlocality and we can simply write C{t) = 
2[r|a|-\/l — a?\qA{t)qB{t)\ — (1 — ''’)/4]- Erom the first line 
of Eq. (Ill with B{t) = Bi{t), it is promptly found that the 
second term under square root is equal to [C{t) -I- (1 — r)/2]^ 
so that the result of Eq. 0 is proven. □ 


IV. EVOLUTION UNDER AN ENVIRONMENT WITH 
OHMIC CLASS SPECTRUM 

We now apply the previous results in the dynamics of en¬ 
tanglement between two initially correlated but noninteracting 
qubits A and B, each undergoing local pure-dephasing due 
to the coupling with a zero-temperature bosonic environment 
having a Ohmic-like spectrum, whose characteristics do not 
depend on the qubit ESll . In this case, the explicit form of the 
Hamiltonian of Eq. 0 ruling the dynamics of each qubit is 
(h = i) 


H = woCTz -I- ^ ujka\ak + cr^ ^ [guak + ^ (13) 

k k 

where ojq is the qubit frequency, Wfc the frequencies of the 
reservoir modes, the Pauli operator along the z-direction, 
ak the bosonic annihilation operators, a\. the bosonic cre¬ 
ation operators and the coupling constants between the 
qubit and each reservoir mode. It is readily seen that the 
above Hamiltonian of Eq. ( [T3] l is obtained by the general pure- 
dephasing Hamiltonian of Eq. 0 by substituting fl = 2 a;o, 

Hr = Y.k^ka\a-k and X = 2 Y,k {gkO-k + ■ In the 

continuum limit one has Q / dujJ{uj)S{uJk—uj), 

where J(w) is the reservoir spectral density which in this case 
is given by 


We stress that this theorem is valid for any local pure- 
dephasing qubit-environment interaction, starting from an ini¬ 
tial EWE state with a generic value of a 7 ^ 0,1. Eor instance, 
when r = 1 Eq. reduces to B{t) = 2-^/1 +C(t)^: the lat¬ 
ter relation, which is known to occur when the system is in a 
pure state 12^ or in a Bell-diagonal state EH, is here found 
to hold during the system evolution for more general states. 

As a consequence of the previous theorem, we get the fol¬ 
lowing 


J{u^) = (14) 

UJc 

with oJc denoting the cut-off reservoir frequency (tOc is as¬ 
sumed equal for the two independent environments). This 
qubit-environment system is exactly solvable in the case of 
zero temperature EOl |2^ . In particular, the single-qubit de¬ 
cay characteristic function q{t) = poi(f)/poi( 0 ) at zero tem¬ 
perature is given by q{t) = with dephasing factor 


Lemma 1. The threshold value Cth of the concurrence such 
that, when C{t) > Cth, the evolved two-qubit state p{t) ex¬ 
hibits nonlocality at time t under any pure-dephasing evolu¬ 
tion and for EWL initial states is 

Cth = Vl - r 2 - (1 - r)/ 2 . ( 12 ) 

Proof The demonstration immediately follows from Eq. 0 
by solving the tight inequality B{f) > 2 with respect to con¬ 
currence C(f). □ 

We name Cth nonlocality threshold. Therefore, for initial 
EWL states evolving under any pure-dephasing interaction, 
the system exhibits nonlocality at time t provided that the con¬ 
currence C (t) is larger than a threshold value Cth depending 
only on the system initial purity. The threshold is a decreasing 
function of the purity and for r = 1 it is Cth = 0 : this implies 
that if the two qubits are initially in a pure entangled Bell-like 
state, their correlations remain nonlocal until any amount of 
nonzero entanglement is present. In the following we exploit 
these results in a specific system where they are particularly 
relevant. 


A(f) = 2 f (15) 

Jo 

where 7 (f) is the dephasing rate of the single-qubit evolution 

El 


7(f) = Wc 1 -f (Wcf) 


-s/2 

r[s] sin [s arctan(a;cf)], 


(16) 


with r[x] being the Euler Gamma function. The parameter s 
of the spectrum of Eq. (14 1 governs the character of the dy¬ 
namics: in particular, for 0 < s < 1 the dynamics is Marko¬ 
vian (memoryless) while for s > 1 it is non-Markovian (with 
memory), as indicated by the time behavior of the single-qubit 
coherence q{t) = Il20ll28]l . In fact, for 0 < s < 1 the 

coherence vanishes asymptotically while for s > 1 it reaches a 
stationary value which is not consistent with a simple Marko¬ 
vian approximation. Moreover, for s > 2 an information 
backflow from the environment to the qubit occurs ll 20 l . 

If the two-qubit are initially prepared in an EWL state of 
Eq. (j^, then the concurrence evolves as in Eq. ( [T0| by substi¬ 
tuting the single-qubit decay characteristic functions qs(f) = 
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g-As(t) where each As(t) is defined in Eq. (15 1 . 

In general, the single-qubit decay function q{t) can be differ¬ 
ent for the two qubits, qA{t), qsit) depending on the values 
of the parameter sa, sb, respectively. In the following we 
shall study the case of both qubits locally subject to identical 
noise (that is, qA{t) = qsit) = q{t)) and the situation where 
a qubit, for instance B, is isolated so that qsit) = 1 at any 
time. 


A first analysis of the entanglement dynamics for some val¬ 
ues of the parameters is displayed in Fig. In particular, for 
both qubits under dephasing we observe that for purity param¬ 
eter r = 0.9 and weights a = l/>/2 (left panel), the entangle¬ 
ment disappears at a finite time, with the faster decay occur¬ 
ring for the larger values of the ohmicity parameter s which 
indicate stronger single-qubit non-Markovianity. Therefore, 
under these initial conditions non-Markovian features do not 
help maintaining entanglement but are detrimental. On the 
other hand, for larger values of the purity of the initial state, 
stationary nonzero entanglement is found, as displayed in the 
right panel of Fig. [^for r = 1. We point out how the case 
where a qubit is isolated is advantageous for the entanglement 
preservation, as expected. Notice also that, for r = 1, the non¬ 
locality threshold for the concurrence obtained from Eq. 0 
is zero, so that this stationary entanglement ensures the pres¬ 
ence of nonlocal correlations. 

The emergence of this stationary entanglement, which did 
not surface in a previous analysis of the same model due to 
the initial conditions there considered 1 ^ , strictly stems from 
the trapping of the single-qubit coherence ll20l . In fact, an 
analytic expression of the dephasing factor of Eq. 0 can be 
obtained by solving the integral which gives 


A(t) = /1 _ _|_ x) ^ [r sin (s arctan r) 

s — 1 1 

-I-cos (s arctanr])} , (17) 

where t = LOct is the dimensionless time in units of uj~^. It 
is then straightforward to see that, in the limit of r to infinity, 
one has 


[ -foo, 0 < s < 1 

A~(s) = lim A(r) = 

i 2r(s)/(s-l), s>l 

(18) 

Therefore, the single-qubit coherence q{t) = decays 

to zero for values of s < 1, while for any value of s > 1 it 
achieves a stationary value q°°{s) = exp{—A°°(s)} depend¬ 
ing on s. A plot of the asymptotic dephasing factor A'^(s) as 
a function of s is displayed in Fig.|^ It is seen that A°°{s) is 
a convex function exhibiting its minimum value at s ~ 2.46 
where A°°(s) ~ 1.77. The value s is thus the value of the 
ohmicity parameter s which maximizes the single-qubit co¬ 
herence Il20l . 

However, the fact that for any s > 1 a stationary coherence 
occurs does not entail that there is correspondingly a nonzero 
stationary entanglement. In fact, as we have already seen in 
Fig. 0 the occurrence of entanglement trapping strongly de¬ 
pends on the value of the purity r of the initial EWE state. The 
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FIG. 1: Dynamics of the two-qubit concurrence starting from Werner states 
(a = b = l/\/2), for different parameter conditions. Upper panel: The 
purity parameter is r = 0.9 for both the three curves and sa = sb = 1-5 
(solid blue line), sa = sb =3 (dashed orange line), sa = 1-5, sb = 
3 (green dot-dashed line). The gray dotted line represents the nonlocality 
threshold for the concuiTence (Cth = 0.386). Lower panel: The purity 
parameter is r = 1 and sa = sb = 2.5 (solid blue line), sa = 2.5 while 
qubit B is isolated (dashed orange line). In this case the nonlocality threshold 
is zero and entanglement trapping is found. 



FIG. 2: Asymptotic value, A°°, of the decay function A(£) as a function 
of the spectrum parameter s. The minimum value of is attained for 
s ~ 2.46 which gives A°°(s) ~ 1.77. 
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analytic expression of the asymptotic concurrence is immedi¬ 
ately obtained by Eq. as 

C“(s) = 2max{0,ray/r^e-”^”(")-(l-r)/4}, (19) 

where A°°{s) is given in Eq. ( p^ and n = 1, 2 for the case 
when one qubit is isolated or both qubits are subject to noise, 
respectively. As seen from the above expression, the value of 
s which maximizes the stationary single-qubit coherence also 
maximizes the stationary concurrence. Eor s < 1 one always 
get C°°{s) = 0. The values of s > 1 such that C°°{s) > 0 
for given a, r can be easily obtained using Eq. ( [T9| ) and are 
determined by the inequality 

A~(s) = ^ \n[AraVl - a7(l - r)]. ( 20 ) 

It is immediately seen that for r = 1 and a ^ 0,1 the 
right-hand side of the above inequality diverges to infinity and 
therefore the latter is satisfied for any value of s > 1. Two 
plots of C°°{s) as a function of s for a = 1 /v^ and two val¬ 
ues of the purity parameter r = 1, 0.99 are shown in Eig[^ 
Eor pure initial states such that r = 1 (left panel) the nonlo¬ 
cality threshold is zero and C°° (s) > 0 for any value of s > 1; 
for r = 0.99 (right panel) the stationary value of concurrence 
is larger than zero for a finite range of s, but only for the case 
when a qubit is isolated this concurrence is able to overcome 
the nonlocality threshold Cth = 0.136 for a smaller range of s 
around s. 

So far, the results obtained tell us that, for s > 1, larger 
values of the initial purity r are more effective in order to get 
stationary entanglement which enables the presence of nonlo¬ 
cal correlations. We can make this qualitative behavior quan¬ 
titative by finding the purity lower bound r* of the initial 
EWE states of Eq. ([^ such that, once assigned the value of 
the weight a and of the ohmicity parameter s, the stationary 
concurrence is above the nonlocality threshold for any r > r*. 
This finding will permit us to know to which degree of purity 
the initial EWE entangled state must be prepared in order to 
guarantee that the two qubits will preserve nonlocality during 
the evolution. Exploiting the nonlocality threshold of Eq. ( [T^ 
together with the stationary concurrence of Eq. we can 
solve the inequality C°° > Cth for s > 1 with respect to r in 
order to obtain the nonlocality lower bound of purity 


1 -f 4 (1 — a^) a? exp{—2nA°°(s)} ’ 


( 21 ) 


where again A°°(s) is given in Eq. (18 1 and n = 1, 2 corre¬ 
sponds to the case when one qubit is isolated or both qubits are 
subject to noise, respectively. Erom this equation, it is easily 
found that when A°° (s) —> -|-c» one has r* = 1, which means 
that for s < 1 the initial entangled state must be pure in order 
to have nonlocal correlations preserved during the dynamics. 
In Eig.j^the plots of r* for n = 1, 2 are displayed as a function 
of s for a = 1 /(left panel) and as a function of the weight 
for s — s (right panel). The plots show the sensitivity of 
the nonlocality purity lower bound r* to the system configu¬ 
ration, with very high values of initial purity of the EWE state 
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FIG. 3: Asymptotic values of two-qubit concurrence starting from Werner 
states {a = b = ll\/2) as a function of the ohmicity parameter s. In both 
panels the lower curve represents the concurrence for the case of both qubits 
subject to noise, while the upper curve is the concurrence for the case where a 
qubit is isolated. Upper panel: The purity parameter is r = 1 and the nonlo¬ 
cality threshold is zero. Lower panel: The purity parameter is r = 0.99 and 
the gray dotted line represents the nonlocality threshold for the concurrence 
(Cth = 0.136). The upper curve, that is the stationary concurrence when 
a qubit is isolated, is larger than the nonlocality threshold within the range 
2.06 < s < 2.94. 


required if both qubits are subject to noise. As expected, from 
Eq.® it is promptly found that the minimum value of r* is 
obtained for s = s and a = 1 /^/2 which are, respectively, the 
value of s minimizing the stationary dephasing factor and the 
value of weight a maximizing the entanglement of the pure 
part of the initial EWE state of Eq. (§. This result is con¬ 
firmed by comparing the left and right panels of Eig. In 
particular, the minimum values of the nonlocality lower bound 
r* for n = 1, 2 are, respectively, = 1) ~ 0.9858 and 

rjjjin(n = 2) ~ 0.9996. We point out that, despite these lower 
bounds of initial purities are very high, they are achievable for 
instance in all-optical setups with entangled polarized photons 
where the open quantum dynamics can be simulated by suit¬ 
able quartz plates ifTTl [T3l 1^ . in solid-state qubits (nuclear 
spins) in diamond and in principle even in composite 
atom-cavity systems by purity swapping ll32l . Values of ini¬ 
tial state purity smaller than the above lower bounds will 
however lead to an entanglement evolution above the nonlo¬ 
cality threshold for long enough dimensionless times Wcf, as 
can be deducted from Eig. 1. 
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FIG. 4: Lower bound r* of the purity parameter r, such that for any 
r > r* the asymptotic concurrence is larger than the nonlocality threshold 
(C°° > Cth)- Upper panel: r* as a function of the ohmicity parameter s for 
\a\ = l/\/2. Lower panel: r* as a function of the weight for ohmicity 
parameter s = s ^ 2.46. The upper dashed curve is for both qubits subject 
to noise while the lower solid one is for the case where a qubit is isolated. 


V. CONCLUSION 

In this paper we have investigated the possibility to main¬ 
tain, during a pure-dephasing evolution, an amount entan¬ 
glement capable to ensure the presence of nonlocal correla¬ 
tions within a two-qubit system. To this purpose, local pure- 
dephasing environments characterized by a Ohmic-like spec¬ 
tral density have been taken into account. We have reviewed 
a known closed formula which dynamically relates the Bell 
function B, witnessing nonlocality, and the concurrence C, 
quantifying entanglement, for arbitrary local pure-dephasing 
evolutions im, giving a formal demonstration of it. This rela¬ 
tion, which permits to find a threshold Cth for the concurrence 
such that for any value of C > Cth the system state exhibits 
nonlocality, has been then exploited within the model with 


Ohmic-like spectrum. 

We have shown that stationary entanglement is achievable 
for a given degree of ohmicity s and for a suitable purity r of 
the initial state. In particular, we have found the analytical ex¬ 
pression of the lower bound of initial purity r* such that, for 
a given s, the stationary entanglement guarantees nonlocality 
for any value of purity of the initial state r > r*. We re¬ 
mark that, although the most performing condition for entan¬ 
glement trapping happens for a value of s, namely s ~ 2.46, 
which enables a non-Markovian regime with the occurrence 
of local information backflows from the environment to the 
system ||20l, the achievement of stationary entanglement is 
simply related to the existence of super-Ohmic environments 
(s > 1) for which the Markovian approximation fails and non- 
Markovian effects arise even without information backflows 
EOl . We also point out that stationary entanglement was ob¬ 
tained for two noninteracting qubits locally subject to ampli¬ 
tude damping channels each with a super-Ohmic spectral den¬ 
sity ll3?l . where the mechanism leading to this phenomenon is 
linked to the presence of a bound state between each qubit 
with its local dissipative bosonic environment. Another re¬ 
cent cavity-based architecture has been reported which is able 
in principle to give stationary entanglement in a dissipative 
compound environment configuration ll?4ll . In the system con¬ 
sidered in this paper, instead, the pure-dephasing interaction 
is non-dissipative and no qubit-environment bound state can 
be created. The origin of the stationary entanglement, activat¬ 
ing nonlocality, is here basically due to the manifestation of 
coherence trapping for the localized qubits. 

The results presented in this work highlight the importance 
of knowing closed relations between different kinds of quan¬ 
tum correlations, as entanglement and nonlocality. Indeed 
this knowledge allows for finding preparation conditions of 
the initial states, under opportune environmental characteris¬ 
tics, which ensure that a given kind of quantum correlations is 
preserved in the long-time limit, for instance nonlocal corre¬ 
lations, just by maintaining another type of correlations, e.g. 
quantum entanglement, above a certain amount. Such find¬ 
ings can count not only for fundamental studies but also for 
practical purposes in quantum communication and informa¬ 
tion processing. 
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